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Abstract 



Several authors have shown that dissipative thermal cycle models based on Finite-Time Ther- 
modynamics exhibit loop-shaped curves of power output versus efficiency, such as it occurs with 
actual dissipative thermal engines. Within the context of First-Order Irreversible Thermodynam- 
ics (FOIT), in this work we show that for an energy converter consisting of two coupled fluxes it is 
also possible to find loop-shaped curves of both power output and the so-called ecological function 
against efficiency. In a previous work Stucki [J.W. Stucki, Eur. J. Biochem. 109, 269 (1980)] used 
a FOIT-approach to describe the modes of thermodynamic performance of oxidative phosphory- 
lation involved in ATP-synthesis within mithochondrias. In that work the author did not use 
the mentioned loop-shaped curves and he proposed that oxidative phosphorylation operates in a 
steady state simultaneously at minimum entropy production and maximum efficiency, by means of 
a conductance matching condition between extreme states of zero and infinite conductances respec- 
tively. In the present work we show that all Stucki's results about the oxidative phosphorylation 
energetics can be obtained without the so-called conductance matching condition. On the other 
hand, we also show that the minimum entropy production state implies both null power output 
and efficiency and therefore this state is not fulfilled by the oxidative phosphorylation performance. 
Our results suggest that actual efficiency values of oxidative phosphorylation performance are bet- 
ter described by a mode of operation consisting in the simultaneous maximization of the so-called 
ecological function and the efficiency. 

PACS numbers: 05.70.Ln Nonequilibrium and irreversible thermodynamics; 84.60.Bk Performance charac- 
teristics of energy conversion systems, figure of merit; 87.16.-b Subcellular structure and processes 
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I. INTRODUCTION 



It is well known that in actual dissipative heat engines, the experimental plots of power 
output against thermal efficiency are loop-shaped curves, where both the maximum power 
and maximum efficiency points do not coincide and their separation can be managed by 
some phenomenological parameters which depend on the engine's design and the materials 
employed in its construction 

a a. 

In these engines is common to find irreversibilities 
(losses) due to friction and irreversible heat fluxes. These losses must be taken into account 
in the models elaborated to describe their general performance. In the thermal engine 
modeling this kind of losses are usually considered in a separate manner. However, one 
can couple the dissipative processes, in such a way that loop-shaped curves are recovered. 
Several authors j^, have proposed thermal engine models which reproduce the loop- 
shaped curves observed in actual engines. This is accomplished by means of characteristic 
functions (as power and efficiency) depending on the thermal reservoir temperatures, and 
other parameters as compression ratios, and thermal conductances for example. Through 
these quantities one can control the distance between the maxima points of both power 
output and efficiency. 

In the present work (Section [IT]), we show that in the case of a linear energy converter 
consisting of two coupled fluxes described by FOIT, in which one spontaneous flow manages 
a nonspontaneous one, loop-shaped curves (LSC) similar to those appearing in irreversible 
thermal engine models can be also obtained. These LSC for our FOIT-model are obtained 
for both power output (P) versus efficiency (77) and for ecological function (E) versus 77, 
with ecological function defined as the power output minus the dissipation function (Ta, see 
subsection IIIBj) [4] . In both loop-shaped curves it appears a force ratio given by Q , 



x — /— !i— (1) 
V L22 X2 

(with Lij being the FOIT-phenomenological coefficients), which measures a direct relation- 
ship between the two forces X\ and X2 involved in the coupled flows. On the other hand, 
we also use a coupling parameter q given by, 

V 

This parameter gives us a measure of the coupling of spontaneous and nonspontaneous 
fluxes 0]. In section HIl we also show that by means of the LSC properties we can study 



several performance modes of the energy converter similar to those used in Finite Time 
Thermodynamics (FTT) ?], S ^, ID], but with further results which consist in a set of func- 



tions describing how the maxima points can move one respect to the other. As an example 
of our previous results, we study the ATP production occurring within mitochondria by 
means of respiration. A FOIT-approach to this problem was previously published by Stucki 



The Stucki's approach was mainly based on using two fluxes and two forces subject 



to several optimization criteria proposed by this author. The Stucki's results regarding the 
economic degrees of coupling of oxidative phosphorylation arise from the assumption that 
mitochondria is in a steady state corresponding to minimum entropy production which si- 
multaneously corresponds to a maximum efficiency state. This last situation stems from 
the so-called conductance-matching condition (CMC), which is obtained by means of the 
inclusion of a third term in the expression for the entropy production, which corresponds to 
an attached cellular load under steady state conditions. The phenomenological coefficient of 
the third term mentioned can be fitted to obtain that the minima values of entropy produc- 
tion coincide with the optimum values of efficiency. Stucki uses several objective functions 
to model different mitochondria modes of operation but only by using the former two forces 
and two fluxes system and by assuming the implicit holding of the conductance-matching 
condition. In the present paper by means of the coupling of only two fluxes and two forces 
some of the main Stucki results are recovered without the inclusion of the third term in en- 
tropy production corresponding to the attached cellular load. This means that the attached 
load is not necessary to describe the energetics of a linear energy converter consisting in a 
pair of coupled processes where one drives the other. In summary, in this article we show 
that the minimum entropy production steady state regime is not equivalent to the optimum 
efficiency steady state regime. This assertion arises from the fact that the Stucki results 
are here obtained without the usage of the attached cellular load. If the Stucki's third term 
in the entropy production is considered, then the energetics formalism should be rewritten 
in terms of a 3 x 3 matrix of phenomenological coefficients, with its consequent changes in 
quantities such as the efficiency and the power output. This paper is organized as follows: 
In section [TT] we discussed on the energetics of a 2 x 2 system; In section IHIt we applied 
some of the results of the previous section on the ATP-production under a linear approach. 
Finally in section IIVI we present some concluding remarks. 
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FIG. 1: Efficiency (77) vs. force ratio (x) for two different fixed coupling parameters. XMr] corre- 
sponds to the maximum of rj for a fixed q. 

II. ON THE ENERGETICS OF A TWO COUPLED FLUXES SYSTEM 



In this section we study the following performance regimes: the minimum dissipation 
function (MDF) Jl2l | , th e maximum power output (MPO) 13j and the maximum ecological 
function (MEF) |l4J. Il5|. By means of these criteria we reproduce some previous results 11] 
and we show some new ones. In particular, we will show that a linear isothermal-isobaric 
engine working in the MEF-regime can reach an efficiency until tjmef = 0.75. 



A. Constitutive equations 

Let J\ and J 2 be two coupled generalized fluxes (being J\ the driven flux and J 2 the 
driver flux); X\ and X 2 are the conjugate generalized potentials associated to the fluxes. 
For the linear case, fluxes and potentials are after Onsager, given by, 



Ji = v^iT (VlTiXi + q^L^X 2 ) (3) 
J 2 = v/Z^ (qy/L^X-i, + y/L^X^j , (4) 
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with L12 = L21, the symmetry Onsager relation between crossed coefficients. In these 
equations we use the coupling coefficient defined by Eq. (j2J) {f|. Thus, in the limit case 
q — > 0, each flux is proportional to its proper conjugate potential through its direct phe- 
nomeno logical coefficient, that is, the crossed effects vanish, and therefore the fluxes become 
independent. When q — > 1, the fluxes tend to a mechanistic stoichiometry fixed relationship 
independently of the potential magnitudes [111 ]. 

On the other hand, it is convenient to define a parameter describing the cross effect 
between both potentials. Taking X 2 > as the associated potential to the driver flux, we 
define the parameter x (see Eq. (DO)) measuring the fraction of X\ < appearing due to the 
presence of the flux J 2 . This parameter is a quantity with values in the interval [—1,0]. To 
complete the set of constitutive equations, we define the efficiency of the thermodynamic 
process as follows, 

energy output 

V = : ~, 

n energy input 

and following Caplan and Essig |6j, in terms of Onsager relations we have, 
^ _ TJ\X\ _ V L u (V LuXi + qy r L^X 2 ) X\ 

TJ 2 X 2 \J L22 (<?V L\\X\ + \J L22X2) X2 

From Eq. (JTJ), we get an expression for 77 in terms of x and q as follows, 

I \ {x + q)x 

r]{x ) q) = — . 6 

qx + 1 

A plot of 7] versus x for a fixed q is depicted in Figure [TJ This graph is a convex curve 
with only a maximum point. That is, there exists a relation between the input and output 
energetic fluxes which maximize the efficiency for some given Onsager coefficients. We take 
now a linear isothermal-isobaric conversion process and build its characteristic functions 
in terms of q and x for determining the performance conditions according to a certain 
objective function maximization. First, we calculate the following characteristic functions: 
the dissipation function, the power output and the ecological function. 



B. Characteristic functions 



Dissipation function- Within FOIT framework, for a 2 x 2 system of fluxes and forces 
the entropy production is given by 0, If]], 



cr = J1X1 + J 2 X 2 , 



(7) 




FIG. 2: Characteristic functions normalized respect to the fixed quantity TL 22 X 2 versus the force 
ratio x. Dissipation function $ (x), power output P (x) and Ecological function E (x), all of them 
for fixed q = 0.97. <& (x) reaches its minimum value at xmdf = —Q, P (x) reaches its maximum 
value at xmpo = — § an d E (x) reaches its maximum value at xmef = 

and following Tribus [5|, the dissipation function for an isothermal system can be expressed 
as $ = To. By the substitution of Eqs. ([3]) and (j3J) into Eq. ([7]), in terms of the parameters 
x and q, $ becomes, 

$ (x, q) = (x 2 + 2xq + l) TL 22 X 2 2 . (8) 

This function has a minimum value at the point xmdf = —q for q and TL22X 2 fixed 
(see Figure [2]), this means that the fraction of X\ due to the presence of X 2 must be 
Ximdf = ~ (ftl) "^ 2 > ^° obtain a minimum dissipation steady state. 

Power output- For isothermal processes of two coupled fluxes the power output is given 




P (x, q) = -TJ X X X = -x (x + q) TL 22 X 2 2 . (9) 

This equation corresponds to a convex curve (Figure EJ), with a maximum value at xmpo — 
— |, this condition implies that the fraction of the driven potential due to driver one in a 
maximum power output regime must be X 1MPS = — ( ^p- ) X 2 . 
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The Ecological function - Defining the ecological function as E = P — $, by means of 
Eqs. © and ®, we get 

E(x,q) = - (2x 2 + 3xq + l)TL 22 Xl (10) 

This equation also corresponds to a convex curve with only a maximum point (Figure 
at xmef — — if> this means that in the maximum ecological regime the relation between 
the driver and driven flows must be X\ MEF = — Ojj^f) X 2 . In this regime the conversion 
process undergoes a pathway accomplishing a good compromise between power output and 
dissipated energy, that is, with a small decrement in P we get a great decrement in Ta 



14 



15|- 



C. Loop shaped plots 



Here, we find LSC by using FOIT-equations in an analogous way as it occurs in FTT- 
models and actual thermal engines. The functions 77 (Eq. [6]), $ (Eq. |SJ), P (Eq. [9} and E 
(Eq. [TO]) depend on two parameters, and three of them (77, P and E) are convex functions 
with respect to x. From the plots corresponding to r\ (Figured]) and P (Figure [2]), we observe 
they have two zeros: when ^JL\\X\ = q\/L 22 X 2 , corresponding to a first order steady state 



n 



[18| . and when y/LuXi << q\J L22X2 corresponding to a totally irreversible energy transfer, 
in both cases 77 = 0. We can transform P, E and $ as functions of 77 {2, ljj. By means of 



Eq. (Q, we first get x (77, q) as 



x(r),q) 



q(l+rj)± \ q 2 (l +77) -4r/ 



(11) 



Here, it is necessary to consider the two solutions of Eq. ffTT]) , because each one represents 
a branch of the plot x versus 77 for fixed q. By the substitution of Eq. ([IT]) into Eq. flS]), we 

(I-77) {2-q[q{l + n)±R]) rrT 2 

z 1 -^22^-2' 



$ (77, q) 



with i? = y q 2 (1 + 77) — 4t/. For the power output (Eq. [H]) we obtain 

pM ^(2- glq ^ + ,)±R ]) TL22Xl 

and in the same way, for the ecological function (Eq. [TO]) , we have 



E (77, q) 



(277 -l)(2- g [g(l + 77)± J R]) 



TL22X2 ■ 



(12) 



(13) 



(14) 
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FIG. 3: Dissipation function (<£) versus efficiency {rj) for 5 = 1, 0.99, 0.98, 0.97, 0.96 respectively. 
r\Mr\ is the point which corresponds to the maximum efficiency for each case (here the case q = 0.96 
is marked). 

When we plot these functions against r] (Figures El H] and El both branches) for q e [q m im 1]) 
we observe that $ has a monotonically decreasing behavior, but P and E describe loop- 
shaped curves with some interesting points: the maxima P and E points (tjmpo an d Vmef) 
and the maximum-77 points (r]Mrj, see Figures H] and [5]) . Plots in Figures H] and [5] are similar 



to those obtained in 



20| for a two-reservoir system with several irreversibilities. We will 



study these conspicuous points in next subsection. 



D. Performance modes and the coupling parameter 



1. Dissipation versus efficiency 

In FigureEl we observe that while q decreases (i.e. the quality of the coupling diminishes), 
the minimum value of $ augments and t]m v (see Eq. (16) below) also decreases. This result 
may be the analogous of that occurring when the heat flux through the body of a motor 
augments preventing that a spontaneous flux can be used in managing a nonspontaneous 
one. In Figure El we see that the minimum dissipation function occurs at t)mdf — 0, where 
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FIG. 4: Power output (P) versus efficiency (77) for the same q's as in Figure 3. t]mpo an d f] s are 
the points corresponding to the maximum power output and the semisum efficiency given by Eq. 
$22]), respectively (here the case q = 0.99 is marked). 

<3> only depends on q as 

®mdf (vmdf) = (1 - q 2 ) TL 22 Xl (15) 

n 

while the power output vanishes Pmdf{vmdf) = [21] . Evidently, when q — 1, we re- 
cover the thermodynamic equilibrium state, where all of the flows vanish, that is, to reach 
&mdf {.Vmdf) = all interactions with the environment must be reversible processes. For 
any other value of g, $mdf (vmdf) ^ 0. 

Another point of interest is that where the efficiency reaches its maximum value, which 
is found by means of d^r)\^ = 0, and lead us to 

q 2 

VMr, (q) = 7 ra , (16) 

(i + v /w) 

with a monotonically decreasing behavior while the irreversibilities increase (q — > 0, see 
Figure M )■ By the substitution of Eq. flTB]) into Eqs. (TT2]). (jT5l) and (jlijl. we obtain 
only functions of q such that $ is monotonically decreasing with q, while P and Z£ are 
convex functions with only a maximum point at qupo — \/2 (\/2 — l) ~ 0.910 and qMEF 
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FIG. 5: Ecological function (E) versus efficiency (77) for the same q's as in Figure 3. tjmef is the 
point which corresponds to the maximum ecological function (here the case q = 0.99 is marked). 



I (\/3 — l) ~ 0.988, respectively (see Figured). The first value was found by Stucki (see 
eq. 52 of [11 1) by maximazing the objective function given by its equation 48 11] (that 



is, the power output) subject to the CMC obtained from his assumption that the attached 
cellular load must be included in the entropy production (the mentioned third term, see 
eqs. 34, 35 and 36 of [ill]). Nevertheless, here we show that this value (qMPO ~ 0.910) 
can be obtained without considering any load, that is, by only optimizing the 2x2 system 
under maximum power output in the maximum efficiency-steady state. Then the minimum 
entropy production steady state is not equivalent to a maximum efficiency steady state (see 
Figure [Tj). The second value (qmef ~ 0.988) has not an equivalent in the Stucki's treatment 
(see below and Figur^T]). 



2. Power output versus efficiency 

The graph of P (rj, q) versus 77 (see Figure HJ) shows that for q = 1 a parabola is obtained. 
However, when we diminish the quality of coupling, loop-shaped curves are obtained, in 
which the relative position between the maxima points of power output and efficiency re- 
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FIG. 6: Comparison between the efficiency of the Mr/-steady state (Eq. (|16|) ). the efficiency of the 
MEF-steady state (Eq. (fT9]) l and the semisum efficiency rj s (Ec. ([22j) . dashed curve). t\mef (q = 1) = 
0.75 is the greatest efficiency of the MEF-regime. 



spectively also depend on the coupling parameter q. Besides, the parabola is the boundary 
mark of the loop-shaped curves, that is, its maximum point gives the ideal power output 
when q = 1, which is P [t]mpo = 1)] = 0.25 x TL 22 X 2 . This result coincides with that 



of ref. 



19] for the muscle contraction problem. 



The efficiency which maximizes the power output is obtained by means of d^Prjl MpQ = 0, 



giving 



1 q 2 
Vmpo — 77; 



(17) 



2 2 - q 2 ' 

Here, we observe that only for q = 1, 7] MPO — ~. At this point the process variables 
corresponding to the MPO-regime satisfy the inequality Pmpo = P [vmpo (<?)] < &mpo = 
$ [Vmpo (<?)], that is, 

-2 



q -TL 22 X\ < ( 1 'i, 2 ) TL 22 X 2 2 V ,/ e [0. 11 



(18) 



as it occurs in the MDF-regime with the advantage that in this MPO-regime the power 
output is not zero, as it is the case for the power output in the minimum entropy production 
steady state [211 ] . 
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FIG. 7: Characteristic functions normalized respect to their maximum values in terms of the cou- 
pling parameter q, at the maximum efficiency regime. For the curve <I> [TjMri(q)] we can observe that 
Qmpo ~ 0.910 (where P [riMri(q)] attains its maximum value) does not correspond to a minimum 
entropy production steady state (with q = 1), that is, the conductance-matching condition to reach 
the MPO-value at r/Mrpregime is not necessary. Power output P [r]Mr](q)] and Ecological function 
E [r/Mr](q)} are also depicted. qMEF ~ 0.988 corresponds to the maximum value of E [rjMriiq)]- 

The maximum efficiency point is found by using dpr]\ p = 0, which also leads to Eq. (1161) . 
and therefore the behavior of the process variables is the same as in the previous case. In 
Figure HI we see that the bigger q the smaller P (i] Mri ), that is, when r\ Mr] increases P {j]mt)) 
decreases until the limit case of the parabola [r\ Mr] = 1 and P = 0). In Figure HJ we also see 
that the distance in the r?-axis between the maximum power and the maximum efficiency 
points diminishes while the quality of the coupling diminishes and reaches a null value when 
q = 0. In the following paragraph we will see how the ecological function gives a good 
tradeoff between maximum power and maximum efficiency. 
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3. Ecological function versus efficiency 



The MPO-regime provides a maximum energy output rate, but also produces a great 
energy dissipation taken from the energy input (low efficiency). Nevertheless, there exist 
some phenomena where this is not observed, that is, the dissipation is always smaller than 



the power output. In fact, many natural processes (bio 
a good compromise between P and $ [g, |7J, 



11 



14 



19 



ogic and nonbiologic) work following 



22 



23 



24,125 



26 



28|]. On the 



other hand, as we see in Figure the ecological function (Eq. (Tl4l) ) plotted versus efficiency 
also gives loop-shaped curves. Therefore, there exists an efficiency for which one obtains 
the best compromise between P and $. This point is found by means of d^El = 0, 
which leads to 



Vmef 



4 4 - 3g 2 ' 

This result gives a q interval where the following inequality now is satisfied &mef 
$ [vmef (<?)] < Pmef = P [vmef (<?)], that is, 



(19) 



1 - ^qA TL 22 Xj < ^TL 22 X^q E I ^, 1 



(20) 



This inequality has some implications about the MEF-regime. In the limit q = 1, t]mef — 
0.75, that is, 0.25 more than the MPO-efficiency, depending on the system's design. Other 
result is that for any value of q, Pmef = 0.75 x Pmpo, while the dissipation function in the 
case of MEF-regime suffers a drastic decreasing compared with §mpo- I n fact we have 

1 /16- 15g 2 ~ 



MEF 



MPO- 



(21) 



4 V 4 - 3q 2 

That is, within the g-interval where E > (q G (^-, 1 ), §mef goes from 0.5 x $a/po to 
0.25 x $mpo To find the maximum efficiency we solve dEf]\ E = 0, obtaining the same 
result that in the two previous regimes, that is, the maximum efficiency has the same value 
for all of the performance regimes given by Eq. (TIB"]) . Therefore, we have only two q values 
that maximize both the power output and the ecological function in the M?7-regime, that 
is, qupo ~ 0.910 and qMEF ~ 0.988, respectively (see Figure [7j). This last value is bigger 
than the biggest (/-value (at maximum efficiency) found by Stucki corresponding to the 
maximization of his function J\X\r\ [ll|. qMEF ~ 0.988 was also found without using the 
CMC. 
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Additionally, in the loop-shaped curves of P versus 77, we observe that it is possible to 
find an intermediate point between maximum power and maximum efficiency accomplishing 
a good compromise between these two ways of performance {4]. This point is 

Vs (?) = \ [vmax + Vmpo] 



1 + VI 



1 



22- q 2 



(22) 



If we compare this expression with the MEF-efficiency, tjmef, we obtain the behavior shown 
in Figure that is, t/mef ~ Vs- This means that the 77 values of the MEF-regime are a 
good compromise between power output and maximum efficiency, which is equivalent to a 
low dissipation regime, with the additional fact that this occurs for a realistic q < 1 (for 
example, Stucki reported a q exp ~ 0.95 for liver mitochondria from male rats However, 
some authors {25] have considered the ideal case q = 1 to study some ATP problems. 



III. ON ATP PRODUCTION: A LINEAR APPROACH 

The purpose of this section is an application of the methodology presented in Section 
Oil One of the most important ex amp les of the energy conversion in biology is the aerobic 



ATP-synthesis (see references 22 
synthesis is given by [2 21 ]. 



24 



25 



26 



27 



281 ] ) . The global chemical reaction of ATP 



{C 6 H 12 6 + 6 2 + 6 H 2 0} + [36 ADP + 36 P + ] 
^ {6 C0 2 + 12 H 2 0} + [36 ATP] , 



(23) 



where driver and driven reactions have been indicated with curly and square brackets, re- 
spectively. For the spontaneous reaction we take J 2 X 2 > and for the nonspontaneous one, 
J1X1 < 0. Some experiments suggest this process occurs out but near an equilibrium state 



|22| , thus, we can use the formalism of Section [Til Stucki considered as a reasonable 
idea that in vivo oxidative phosphorylation simultaneously operates at both maximum ef- 
ficiency and minimum entropy production. This situation in Stucki's words is reached in a 
steady state named conductance matching between extreme states of zero and infinite con- 



ductances respectively 



111 ] . However, as it can be observed in Figure H] (which stems from 



the parametric combination of Eqs. (jSJ) and (jUJ) for realistic values of \q\ < 1), there exists a 
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unique point with simultaneous zero values for r\ and P (which corresponds to the minimum 
entropy production steady state, see Figure [3]). Thus, for q < 1 a conductance-matching 
condition in the Stucki's sense is not necessary. Among the four objective functions pro- 



posed in 



11], that given by F$ = JiXit] at maximum efficiency, is in Stucki's words the 



most suitable for the oxidative phosphorylation case. The optimal efficiency obtained with 
this function is rjp = 0.618 arising from a coupling parameter of q e ^ ~ 0.972 (his Eq. (58)). 
This Tjp value is lower than reported efficiencies calculated from actual free energy changes, 
which are larger than their corresponding standard free energy changes for biochemical re- 



actions as Eq.(l23l 



be around 0.736 



31 



n fact, for ATP synthesis under in vivo conditions the efficiency can 
32j | . Remarkably, these values rjp = 0.618 and q e £ ~ 0.972 can be 



also obtained with the 2x2 energy conversion formalim without using the CMC, in the 

t Eq. (USD 



33| in the 



following way: by the substitution of r]M V given by Eq. ([16]) into the power output 
and multiplying again by r\ Mr] we get F$ {j]Mr) (the so-called "efficient power" 
maximum efficiency steady state), which corresponds to a convex curve with a maximum at 
q ~ 0.972 (i.e. q e £\ By substituting this g-value into Eq. (1161) we obtain the same rjp = 0.618 
(without CMC). Thus, we obtained these results for ijMri and gff, and that of subsubsection 
III D II for q p pa 0.910 by using only a 2 x 2 formalism. In fact, all of the Stucki results of his 



Table I [ljj for his four economic degrees of coupling can be obtained by the same procedure 
without recurring to the CMC. 

If we take the efficiency value provided by the FOIT-formalism as a good criterion to 
choice the objective function at which oxidative phosphorylation thermodynamically per- 
forms, then we can propose as objective function the ecological one due to its properties pre- 
viously discussed. The efficient power ( F$ = J\Xirj) used by Stucki [ll| leads to rjp = 0.618. 
However, Nelson and Cox (see Box 13.1 of [3l|], pag. 498) reported that oxidative phospory- 
lation under in vivo conditions can reach efficiencies as high as 0.736 due to actual free 
energy changes are larger than their corresponding standard free energy changes. If one 
observes Figure HI along the LSC between the maximum power output point and the maxi- 
mum efficiency point, one has an infinite number of points corresponding each one of them 
to a particular mode of performance. Among these points one can choice some of them in 
terms of their energetic properties, such as a good compromise between high power output 



and low dissipation [4j, |l7|. One mode of performance that accomplishes this goal is the 
so-called ecological function (see subsubsection HID 31) . If we use the ecological function 
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at maximum efficiency, which is analogous to the "economic" functions of Stucki, we find 
for the MEF-regime, qMEF ~ 0.988 (see subsubsections HTD II and III D 31) . If we substitute 
this value into Eq. (|T5|) (the expression for the maximum efficiency) we get r\ MEF = 0.732, 



which is near estimated actual efficiency values 



31 



32]. To compare the thermodynamic 



performance of the ecological function {(Jmef ~ 0.988) with the efficient power (q^ c m 0.972) 



we use a compromise function C (q) of the type defined in 17|. This function in terms of 



normalized quantities respect to the MPO-regime at optimum efficiency is given by 

c , n = P lVM V (g)j _ g [Vm v (q)] , 2A , 

P [vm v (qMPo)} $ [vm v (qMPo)] 

The function C (q) has a maximum at qc ~ 0.982, and it represents the best compromise 
between high power output and low dissipation. If this function is evaluated in qff, then 
C(qp C ) = 0.479. On the other hand, for qMEF, C (qMEF) = 0.489, that is, a slightly larger 
value than C (q!p C ) . However, in percentage terms qMEF is twice closer to qc ~ 0.982 than q^. 
Thus, the ecological optimization provides a reasonable criterion for the thermodynamical 
performance of oxidative phosphorylation (without CMC), with the advantage of giving a 
high efficiency value within the range of actual values. 



IV. CONCLUDING REMARKS 

In the present paper we have developed a procedure to study an irreversible linear en- 
ergy converter working under several steady-state conditions, by using optimization criteria 



stemming from some non-equilibrium approaches 29] . Such criteria enhance the information 
given by the FOIT-formalim. In this way, one can describe the mutual influence between 
the operating mode of the energy converter (generalized forces and fluxes) and its design 
(phenomenological Onsager coefficients). Thus, we get some insights about the quantitative 
description of the energetics of a 2 x 2 system. 

Our results permit to see that the FOIT-formalism along with FTT-procedures lead to 
loop-shaped curves such as it is observed in actual dissipative thermal engines j^, 3] and some 
biological systems as it is the case of experimental data of efficiency versus power output 
reported by Smith et al [24| for the soleous muscle of mouse. Our thermodynamic approach 
;o the bridge between FOIT and some concepts arising from finite-time thermodynamics 



29| | are somewhat different to that suggested by Verhas and de Vos 



17 



Finally, in this work we have obtained all of the Stucki's results for the optimal 
efficiency and the economic degrees of coupling of oxidative phosphorylation by using only 
a 2 x 2 coupled system of fluxes without resorting the so-called conductance matching 
condition. Therefore, our approach suggests that the minimum entropy production regime 
is not compatible with the optimal efficiency steady state. In fact the MEP-regime leads 
to both zero efficiency and zero power output. On the other hand, we found that the 
so-called maximum ecological regime is suitable for the energetic description of oxidative 
phosphorylation, since this criterion gives a high efficiency value within the range of actual 
efficiencies. In addition, this regime represents a good compromise between high power 
output and low dissipation. In summary, our results indicate that the role played by the 
attached cellular load is not necessary for the energetic description of the two coupled fluxes 
involved in the biochemical reaction of ATP synthesis, in the same way that no particular 
load is necessary to describe the internal energetics properties of a typical power plant. If 
one wishes to take into account the external load, a 3 x 3 formalism is necessary for the 
overall thermodynamic description of the complete system, leading to a new set of energetic 
equations different to those as eqs. (jHJ) and (jHJ). 
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